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ABSTRACT 

We present V and £?-band CCD photometry of 35 RR Lyr stars in the southern 
hemisphere globular cluster NGC 6362. Fourier decomposition of the light curves was 
used to estimate the basic properties of these variables. From the analysis of the RRc 
stars we obtain a mean mass of M = 0.531 ± 0.014 Mq, luminosity logL/L© = 
1.656 ± 0.006, effective temperature T e // = 7429 ± 20, and helium abundance Y = 
0.292 ± 0.002. The mean values of the absolute magnitude, metallicity (on Jurcsik's 
scale) and effective temperature for 14 RRab stars with "regular" light curves are: 
M v = 0.86 ± 0.01, [Fe/H] = -0.93 ± 0.04 and T eff = 6555 ± 25, respectively. From 
the B — V colors, periods and metallicities of the RRab stars we estimate the color 
excess for NGC 6362 to be equal to E{B — V) = 0.08 ± 0.01. Using this value we 
derive the colors of the blue and red edges of the instability strip in NGC 6362 to be 
[B - V) BE = 0.17 and {B - V)§ E = 0.43. From the relations between the Fourier 
coefficients of RRab and RRc stars and their absolute magnitudes we estimate the 
apparent distance modulus to NGC 6362 to be (m — M)v = 14.46 ± 0.10. From the 
mean value of log L/Lq of the RRc stars we obtain 14.59 ± 0.03. 

The U-band light curves of three of the RRc stars exhibit changes in amplitude 
of over 0.1 mag on the time scale of around one week. Near the radial 1st overtone 
frequency we find one or two peaks, which strongly suggests that these variables belong 
to the newly identified group of non-radial pulsating RR Lyr stars. 

Key words: stars: RR Lyr - stars: variables - globular clusters: individual: NGC 
6362 



1 INTRODUCTION 

The Clusters AgeS Experiment (CASE) is a long term 
project with the main goal of determining accurate ages and 
distances of globular clusters by using observations of de- 
tached eclipsing binaries (Paczyriski 1997). NGC 6362 was 
selected as one of the targets of our ongoing photometric sur- 
vey aimed at identification of binaries suitable for detailed 
analysis (Kaluzny et al. 1999; Thompson et al. 1999). As a 
byproduct we obtain time series photometry of RR Lyr stars 
located in the surveyed clusters (Olech et al. 1999a, 1999b, 
Kaluzny et al. 2000) 

NGC 6362 (C1726 - 670) is one of the southern hemi- 
sphere Oosterhoff type I globular clusters. The low central 
concentration of the cluster and its proximity enable the 



measurement of reliable photometry for the majority of its 
stars. The Third Catalogue of Variable Stars in Globular 
Clusters (Sawyer Hogg 1973) contains a list of 33 variables 
in NGC 6362. An updated list of these stars with more pre- 
cise period determinations was published by Clement et al. 
(1995). Three objects from this list (V4, V9 and V32) were 
noted as non variable stars. The nature of VI was listed as 
unknown and the other objects were classified as RR Lyr 
stars. 

Recently Mazur et al. (1999) published results of a CCD 
search for short-period variables in NGC 6362. They discov- 
ered 19 new candidates, five of which were cluster RR Lyr 
stars, four were SX Phe stars and eight were eclipsing bina- 
ries. 

The metallicity of NGC 6362 is moderately high and 
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relatively well established. On the scale of Zinn and West 
(1984), [Fe/H] is equal to —1.08. Recently, Carretta and 
Gratton (1997) obtained [Fe/H] = —0.96 from an analysis 
of echelle spectra of cluster giants. 

Alcaino and Liller (1986) obtain a distance modulus 
(rn - M)v = 14.74 and a reddening E(B - V) = 0.10. 
Recently, Piotto et al. (1999) obtained (m-M) v = 14.67 ± 
0.20 and E(B — V) = 0.06 ±0.03 from an analysis of Hubble 
Space Telescope data. 



2 OBSERVATIONS AND DATA REDUCTION 

Time series photometry of NGC 6362 was obtained dur- 
ing the interval 1999 April 17 - August 22 with the 1.0-m 
Swope telescope at Las Campanas Observatory. We used the 
backside-illuminated SITE3 2048 x 4096 CCD, with a pixel 
scale 0.435 arcsec/pixel. A 2048 x 3150 subsection of the 
CCD was used giving a field of view of 14.8 x 22.8 arcmin 2 . 
Photometry presented in this paper is based on 1088 V- 
band images and 259 B-band images which were collected 
on 42 nights. Exposure times were adjusted on every night 
- depending on the seeing - to assure unsaturated images 
for stars from the horizontal branch of the cluster. For the 
l/-band these ranged from 180 to 300 seconds. 

Instrumental photometry was extracted using DoPHOT 
(Schechter, Mateo and Saha 1993). We used DoPHOT in 
the fixed-position mode, with the stellar positions measured 
on "template" images constructed by averaging the 3 best 
exposures obtained for each filter. The instrumental pho- 
tometry was transformed to the standard BV system using 
observations of standard stars from Landolt (1992). We es- 
timate that systematic errors of the zero points of the final 
BV photometry do not exceed 0.03 mag. 



3 RESULTS 

Our search for pulsating variable stars in NGC 6362 has 
identified 35 RR Lyrae stars. All these variables were pre- 
viously known (Clement et al., 1995, Mazur et al., 1999, 
we use the nomenclature from these papers). The periods 
P, V filter peak-to-peak amplitudes Av, intensity averaged 
< V > and < B > magnitudes and types of these stars are 
presented in Table 1. The phased light curves of these stars 
are shown in Fig. 1. Crowding prevented the classification 
of VI and V32 by Clement et al. and Mazur et al. The light 
curves in Fig. 1 indicate that both are RRab stars. 

We compared our photometry with the observations of 
Mazur et al. (1999). For 18 RR Lyr stars in common, the 
mean difference < V > our — < V > Mazur is equal to 0.001 ± 
0.006. 

In Fig. 2 we present a color-magnitude diagram for RR 
Lyr variables in NGC 6362. Open circles denote the RRc 
stars, and open triangles correspond to the RRab variables. 

The periods of the cluster RRc stars are between 0.2462 
and 0.3815 days with a mean value of 0.295 days. The pe- 
riods of the RRab variables are between 0.4471 and 0.6983 
days with a mean value of 0.538 days. These values are in 
very good agreement with the mean values of periods of RRc 
and RRab stars in NGC 6362 of 0.296 and 0.556 days, re- 
spectively, obtained by Clement et al. (1995). These mean 



Table 1. Elements of the pulsating variables in NGC 6362 



Star 


P [days] 


A v 


< V > 


< B > 


T unr 

1V P L 


VI 





.504800 


1 


.09 


15 


.274 


15. 


.657 


RRab 


V2 





.488973 


1 


.27 


15 


.243 


15. 


.595 


RRab 


V3 





.447110 


1 


,14 


15 


,266 


15. 


.598 


RRab 


V5 





.520730 





.74 


15 


.332 


15. 


.741 


RRab 


V6 





.262703 





.47 


15 


.288 


15. 


.547 


RRc 


V7 





.521590 


1 


,03 


15 


.276 


15. 


.664 


RRab 


V8 





.381480 


0. 


.44 


15 


.087 


15. 


.362 


RRc 


V10 





.265675 





.17 


15 


.296 


15. 


.552 


RRc 


Vll 





.288790 





.50 


15 


.219 


15. 


.475 


RRc 


V12 





.532905 


1 


.02 


15 


.230 


15. 


.581 


RRab 


V13 





.579987 


1 


.00 


15 


.211 


15. 


.573 


RRab 


V14 





.246210 





.34 


15 


.337 


15. 


.588 


RRc 


V15 





.279940 





.43 


15 


.221 


15. 


.540 


RRc 


V16 





.525668 


1 


.08 


15 


.299 


15. 


.706 


RRab 


V17 





.314600 


0. 


,12 


15 


.308 


15. 


.639 


RRc 


V18 





.512876 





.96 


15 


.306 


15. 


.706 


RRab 


V19 





.594506 





.59 


15 


.318 


15. 


.794 


RRab 


V20 





.698330 


0. 


.38 


15 


.206 


15. 


.713 


RRab 


V21 





.281390 


0. 


.49 


15 


.267 


15. 


.558 


RRc 


V22 





.266840 





.51 


15 


.317 


15. 


.582 


RRc 


V23 


0. 


.275110 





.49 


15 


.303 


15. 


.574 


RRc 


V24 





.329360 


0. 


.46 


15 


.109 


15. 


.432 


RRc 


V25 





.455891 


1 


.28 


15 


.407 


15. 


.722 


RRab 


V26 





.602182 


0. 


.61 


15 


.327 


15. 


.789 


RRab 


V27 





.278120 





.51 


15 


.280 


15. 


.545 


RRc 


V28 





.358470 





.11 


15 


.074 


15. 


.406 


RRc 


V29 





.647743 


0. 


,12 


15 


.218 


15. 


.688 


RRab 


V30 





.613372 


0. 


.94 


15 


.137 


15. 


.542 


RRab 


V31 





.600230 





.60 


15 


.294 


15. 


.739 


RRab 


V32 





.497203 


1. 


.11 


15 


.277 


15. 


.659 


RRab 


V33 





.306420 


0. 


.42 


15 


.318 


15. 


.636 


RRc 


V34 





.494370 


0. 


.68 


15 


.326 


15. 


.723 


RRab 


V35 





.290790 





43 


15 


.254 


15. 


.576 


RRc 


V36 





.310090 





.38 


15 


.107 


15. 


.475 


RRc 


V37 





.255062 


0. 


.28 


15 


.273 


15. 


.478 


RRc 



values are also similar to those obtained for the Oosterhoff 
type I clusters NGC 6171 and NGC 6732 with characteristics 
resembling NGC 6362 (Sandage 1993). 

We fitted our l/-band light curves with Fourier sine se- 
ries of the form: 

12 

V = Ap + Aj ■ sin(jarf + fa) (1) 

3=1 

where u) = 2-n/P. To find the values of uj, Aj and fa we 
employed the method developed by Schwarzenberg-Czerny 
(1997) and Schwarzenberg-Czerny and Kaluzny (1998). To 
demonstrate precision of our light curves, in Fig. 3 we plot 
the observed light curve of RRab variable V25 together with 
the Fourier sine series with j = 15. The observational points 
are plotted using grey open circles of diameter 0.06 mag, 
equal to 10 standard deviations. The solid black line cor- 
responds to the fit based on equation (1). Table 2 presents 
the amplitudes Aj and phases fa and their errors. For RRab 
stars with amplitudes of 1.25 mag we can determine values 
of Aj up to j = 14 with relative error smaller than 10%. 
The relative error of Ais is equal to 14%. 

The values of the peak-to-peak amplitudes Av and pe- 
riods P presented in Table 1 are used to plot the period- 
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Table 2. Fourier Elements of the RRab variable V25 from NGC 
6362 



0.4401 


aA\ 
0.0003 


A 2 
0.2118 


0A2 
0.0003 


A 3 
0.1576 


0-A3 
0.0003 


A 4 
0.1010 


crA 4 

0.0003 


A 5 
0.0737 


<tA 5 
0.0003 


A 6 
0.0508 


crAjj 

0.0003 


A 7 
0.0320 


crA 7 

0.0003 


As 
0.0211 


crA 8 
0.003 


A 9 
0.0173 


0-A9 
0.0003 


A w 
0.0137 


ctAio 
0.0003 


An 
0.0097 


crAn 

0.0003 


A12 
0.0059 


0-A12 
0.0003 


Al 3 

0.0035 


<tA 13 
0.0003 


A14 

0.0029 


0-A14 

0.0003 


A 1S 
0.0022 


0A15 

0.0003 


01 
5.169 


0.001 


0.041 


<T(p2 

0.002 


fo 

1.506 


CT 03 

0.002 


04 

3.007 


(T04 

0.003 


05 

4.538 


CT05 

0.005 


06 

6.129 


(T06 

0.007 


07 

1.327 


(707 

0.010 


08 

2.674 


CT08 

0.016 


09 

4.080 


(T0g 

0.019 


010 

5.556 


0-010 
0.024 


011 

0.875 


(7011 

0.034 


012 

2.362 


C012 

0.055 


013 

3.683 


CT013 

0.094 


014 

5.198 


(T014 

0.117 


015 

0.100 


CT015 

0.155 


amplitude (log 10 


P - A v 


) diagram 


shown 


in Fig. 4. 



circles denote RRc stars, filled triangles RRab stars with 
D m < 3 (for definition of D m see section 3.3) and open 
triangles RRab variables with D m > 3. The solid line rep- 
resents a linear fit to RRab variables in M3 (Kaluzny et al. 
1998) another Oosterhoff type I globular cluster. 

Clement and Shelton (1999) have suggested that for 
RRab stars the V amplitude at a given period is not a func- 
tion of metal abundance, but rather only a function of the 
Oosterhoff type. We see in Fig. 4 that the majority of our 
RRab variables lay below the mean line for Oosterhoff type I 
clusters, as defined by the RRab stars in the globular cluster 
M3. A similar result is also found by Kaluzny et al. (2000) for 
the cluster M5 and by Wehlau, Slawson, and Nemec (1999) 
for the cluster NGC 7006. These results suggest that while 
there is a clear correlation of Oosterhoff type with metal- 
licity (see Table 4), within at least the Oosterhoff type I 
clusters there is an additional dependence of V amplitude 
at a given period upon metallicity, with the more metal rich 
stars having lower V amplitudes. 

3.1 Fourier Analysis of the RRc Stars 

Simon and Clement (1993) showed that the Fourier decom- 
position of the light curves of RRc variables is a very useful 
technique for estimating the physical parameters of these 
stars. The equations of Simon and Clement (1993) are: 

log M = 0.52 log Pi - 0.11^1 +0.39 (2) 

log L — 1.04 log Pi - O.O5805i + 2.41 (3) 

logT e // = 3.265-0.3026 log Pi -0.1777 log M+0.2402 log L(4) 

logF = -20.26+4.935 log T e// -0.2638 log M+0.3318 log L(5) 



where M is the mass of the star in solar units, Pi is the 
first overtone pulsation period in days, L is the luminosity 
in solar units, P e // is the effective temperature in Kelvin, 
Y is the relative helium abundance and 03i = 03 — 30*. 
The phases marked by asterisks are obtained from a cosine 
Fourier series (used by Simon and Clement 1993) and differ 
from our phases which were obtained from a sine series (cf. 
equation 1). For 03i we have 03i = <03i + tt. 

In addition, using the formula of Kovacs (1998) we may 
estimate the absolute magnitude My° of RRc stars directly: 

My" = 1.261 - 0.961Pi - 0.004021 - 4.447A4 (6) 

From Eqs. (2)-(6) we compute masses, luminosities, ef- 
fective temperatures, relative helium abundances and ab- 
solute magnitudes of RRc stars in NGC 6362. These are 
presented in Table 3 together with the values of Ao, Ai, <f>2i 
and 03i . The errors presented in Table 3 are calculated from 
the formal errors of the Fourier coefficients using the error 
propagation law. 

We exclude from our sample variables V6, V10 and V37 
due to their irregular light curves and large errors of 03i (see 
section 3.2). The mean values of the mass, luminosity, effec- 
tive temperature and helium abundance for the remaining 
14 RRc stars are 0.531+0.014 M Q , log L/Lq = 1.656+0.006, 
T eff = 7429 ± 20, and Y = 0.292 ± 0.002, respectively. The 
mean log Pi is equal to —0.525 ± 0.014 and the mean 03i to 
3.582 + 0.162. 

By comparing the physical parameters of the RRc stars 
from different clusters Clement and Shelton (1997) revealed 
several correlations. Specifically, an increasing mean mass 
corresponds to an increasing luminosity and a decreasing 
effective temperature and relative helium abundance. It is 
clearly visible from Table 4 taken from Kaluzny et al. (1998, 
2000) that NGC 6362 fits this sequence rather well. 

3.2 Non-radial Pulsations of RRc Stars 

Three of the RRc variables presented in Fig. 1, namely, V37, 
V6 and V10, show clear modulation of their light curves. 
The amplitudes and shapes of the light curves change on 
the time scale of days. This is not the Blazhko effect, which 
is observed mainly in RRab stars as a modulation of the 
amplitude on a time scale of tens of days. The other pos- 
sible explanation is that these stars are bimodal pulsators 
(RRd stars) with pulsation in both the fundamental mode 
and the 1st overtone. To check this possibility we calculated 
the power spectra of these variables. The power spectrum 
of the light curve of V37 is presented in the upper panel 
of Fig. 5. The highest peak is detected at a frequency of 
/1 = 3.920615 cyc/day (Pi = 0.2550062 d ). The light curve 
of this star phased with Pi is shown in the third plot in 
the upper panel. One can see that the light curve is noisy. 
The power spectrum shows a secondary peak at frequency 
/ 2 = 3.984365 (P 2 = 0.250981 d ). We prewhitened our orig- 
inal light curve removing the main period Pi with its har- 
monics. The power spectrum of the prewhitened light curve 
is shown in the middle panel of Fig. 5. There is no trace of 
the main frequency /1 and the highest and only significant 
peak is /2. The third plot in the middle panel presents the 
prewhitened light curve phased with the period P2. Again we 
prewhitened this light curve now removing period P2 with 
its harmonics, and the resulting power spectrum shows no 
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Table 3. Parameters for the RRc Variables in NGC 6362 



Star 




A 




Ai 




021 




031 




M 


1( 


5gL 




Y 






My 


V6 


15 


.304 


0. 


.244 


3. 


166 


6. 


476 





.526 


1 


.613 


7544 


0.305 


0.825 





.753 




±0 


.000 


±0. 


.003 


±0. 


.105 


±0. 


190 


±0 


.025 


±0 


.011 


±18 


±0.006 


±0.014 






V'8 


15 


.096 





.217 


1. 


182 


7. 


778 


0. 


.460 


1 


.706 


7267 


0.282 


0.657 





.520 




±0 


.000 


±0 


.001 


±0 


.084 


±0. 


032 


±0 


.004 


±0 


.002 


±2 


±0.001 


±0.006 






V10 


15 


.311 


0. 


.242 


3. 


151 


6. 


154 





.574 


1 


.637 


7500 


0.295 


0.818 





.693 




±0 


.000 


±0 


.002 


±0. 


.070 


±0. 


136 


±0 


.020 


±0 


.008 


±13 


±0.004 


±0.009 






Vll 


15 


.234 





.254 


3. 


315 


6. 


614 


0. 


.534 


1 


.648 


7454 


0.294 


0.762 





.665 




±0 


.000 


±0 


.001 


±0 


.017 


±0. 


027 


±0 


.004 


±0 


.002 


±2 


±0.001 


±0.005 






V14 


15 


.343 


0. 


.169 


3. 


047 


6. 


090 





.561 


1 


.606 


7577 


0.305 


0.881 





.770 




±0 


.000 


±0 


.000 


±0 


.015 


±0. 


068 


±0 


.010 


±0 


.004 


±6 


±0.002 


±0.001 






V15 


15 


.234 





.219 


2. 


959 


6. 


201 


0. 


.583 


1 


.657 


7448 


0.289 


0.804 





.643 




±0 


.000 


±0 


.000 


±0 


.016 


±0. 


027 


±0 


.004 


±0 


.002 


±2 


±0.001 


±0.005 






V17 


15 


.316 


0. 


.210 


3. 


172 


7. 


267 





.473 


1 


.648 


7425 


0.298 


0.788 





.665 




±0 


.000 


±0 


.001 


±0 


.033 


±0. 


044 


±0 


.005 


±0 


.003 


±4 


±0.001 


±0.005 






V21 


15 


.281 





.253 


3. 


173 


6. 


398 


0. 


.557 


1 


.648 


7461 


0.293 


0.775 





.665 




±0 


.000 


±0 


.000 


±0 


.009 


±0. 


013 


±0 


.002 


±0 


.001 


±1 


±0.000 


±0.000 






V22 


15 


.317 


0. 


.254 


3. 


156 


5. 


972 





.603 


1 


.649 


7478 


0.290 


0.790 





.663 




±0 


.000 


±0 


.000 


±0 


.007 


±0. 


017 


±0 


.003 


±0 


.001 


±1 


±0.000 


±0.000 






V23 


15 


.305 





.253 


3. 


134 


6. 


086 


0. 


.595 


1 


.656 


7456 


0.288 


0.779 





.645 




±0 


.000 


±0 


.000 


±0 


.013 


±0. 


021 


±0 


.003 


±0 


.001 


±1 


±0.001 


±0.001 






V24 


15 


.121 





.231 


3. 


329 


7. 


303 


0. 


.480 


1 


.667 


7379 


0.292 


0.740 





.618 




±0 


.000 


±0 


.001 


±0. 


.054 


±0. 


064 


±0 


.008 


±0 


.004 


±6 


±0.002 


±0.005 






V27 


15 


.296 


0. 


.261 


3. 


091 


6. 


195 





.582 


1 


.655 


7455 


0.289 


0.773 





.648 




±0 


.000 


±0 


.001 


±0 


.014 


±0. 


022 


±0 


.003 


±0 


.001 


±2 


±0.001 


±0.004 






V28 


15 


.085 





.221 


3. 


684 


7. 


562 


0. 


.470 


1 


.690 


7313 


0.286 


0.692 





.560 




±0 


.000 


±0 


.001 


±0. 


.045 


±0. 


.031 


±0 


.004 


±0 


.002 


±2 


±0.001 


±0.005 






V33 


15 


.328 


0. 


.215 


3. 


161 


7. 


022 


0. 


.497 


1 


.651 


7429 


0.296 


0.796 





.658 




±0 


.000 


±0 


.000 


±0. 


.020 


±0. 


029 


±0. 


.004 


±0 


.002 


±2 


±0.001 


±0.001 






V35 


15 


.267 


0. 


.217 


3. 


038 


6. 


525 


0. 


.548 


1 


.656 


7438 


0.291 


0.803 





.645 




±0 


.000 


±0. 


.001 


±0. 


.021 


±0. 


033 


±0. 


.005 


±0 


.002 


±3 


±0.001 


±0.005 






V36 


15 


.113 


0. 


.189 


3. 


178 


7. 


113 


0. 


.488 


1 


.651 


7425 


0.296 


0.788 





.658 




±0 


.000 


±0 


.001 


±0 


.044 


±0. 


058 


±0 


.007 


±0 


.003 


±5 


±0.002 


±0.005 






V37 


15 


.277 


0. 


.138 


2. 


034 


4. 


.061 


0. 


.956 


1 


.740 


7343 


0.252 


0.904 





.435 




±0 


.000 


±0 


.003 


±0 


.199 


±5. 


529 


±1. 


.339 


±0 


.321 


±520 


±0.010 


±0.016 







Table 4. Mean parameters for RRc stars in globular clusters after Kaluzny et al. (1998, 2000). 



Cluster 


Oostcrhoff 


[Fe/H] 


No. of 


mean 


mean 


mean 


mean 




type 




stars 


mass 


logL 


T eff 


Y 


NGC 6171 


I 


-0.68 


6 


0.53 


1.65 


7447 


0.29 


NGC 6362 


1 


-1.08 


14 


0.53 


1.66 


7429 


0.29 


M5 


I 


-1.25 


14 


0.54 


1.69 


7353 


0.28 


M3 


1 


-1.47 


5 


0.59 


1.71 


7315 


0.27 


M9 


11 


-1.72 


1 


0.60 


1.72 


7299 


0.27 


M55 


11 


-1.90 


5 


0.53 


1.75 


7193 


0.27 


NGC 2298 


11 


-1.90 


2 


0.59 


1.75 


7200 


0.26 


M68 


11 


-2.03 


16 


0.70 


1.79 


7145 


0.25 


M15 


11 


-2.28 


6 


0.73 


1.80 


7136 


0.25 



further significant frequencies (see the lower panel of Fig. 
5). We conclude that both frequencies correspond to real 
modes of pulsation. For ordinary double mode pulsators the 
ratio of the 1st overtone period to the fundamental period is 
approximately 0.745. In our case the ratio P2/P1 is equal to 
0.98 and thus it is clear that if Pi corresponds to the radial 
1st overtone pulsations then P2 in V 37 corresponds to a 
non-radial mode. 

Another variable with a "noisy" light curve is V6. The 



power spectrum of its raw light curve is presented in the 
upper panel of Fig. 6. The highest peak has frequency /1 = 
3.806580 cyc/day (Pi = 0.262703 d ). The power spectrum 
reveals the presence of other components. We prewhitened 
the light curve of V6 using Pi, and the power spectrum 
of this prewhitened light curve is presented in the second 
panel of Fig. 6. There are two high peaks with frequencies of 
h = 3.740961 (P 2 = 0.267311 d ) and f 3 = 3.871003 cyc/day 
(P3 = 0.258331 d ) in this plot. We again prewhitened our 
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light curve using the period P2 and its harmonics, and the 
power spectrum of the twice prewhitened light curve is pre- 
sented in the third panel of Fig. 6. The only peak that is left 
is /3. After removing P3 with its harmonics we obtained a 
power spectrum clean of any significant features, presented 
in the lower panel in Fig. 6. 

The frequencies fz and f% lay symmetrically above and 
below fx, with f 3 - fx = 0.0644 and /1 - f 2 = 0.0656. This 
symmetry suggests that P2 and P3 are signatures of the am- 
plitude modulation with a period equal to l/(/i — f 2 ) = 
l/(/3 — fx) ~ 15 days. This resembles the Blazhko ef- 
fect. Although the Blazhko effect is common among RRab 
stars, it is seldomly observed in RRc stars. In particu- 
lar, in a sample of 46 stars with known Blazhko periods 
there are only three RRc stars, and the shortest Blazhko 
period is 10.9 days (Smith 1995). However, the difference 
(fx — f 2 ) — (/3 — fx) = 0.0012 is large compared to the ac- 
curacy of our frequency determinations. Since the separate 
light curves phased with periods Pi , P2 and P3 shown in Fig. 
6 display stable amplitudes, they could correspond to pul- 
sations in three different modes of which two are non-radial, 
similar to the one non-radial mode in V37. 

A similar analysis was performed for variable V10 with 
the results presented in Fig. 7. In this case we also detected 
three periods. The most significant is Pi = 0.265675 d (/1 = 
3.763997) and two other are P 2 = 0.257573 d (f 2 = 3.882395) 
and P 3 = 0.274262 d (/ 3 = 3.646148). Again the differences 
of frequencies f 2 — fx and fx — fz are similar and equal to 
0.1184 and 0.1178, respectively. The corresponding beat pe- 
riod l/(/2 — fx) — l/(/i — h)is ~ 8.5 days. If we assume 
that the amplitude modulation is the Blazhko effect, corre- 
sponding Blazhko period would be the shortest known! For 
V10 the difference (/ 2 - fx) - (fx - /a) is equal to 0.0006 
and, given the stable light curves in Fig. 7, we suggest that 
the periods Pi and P3 are real and correspond to non-radial 
pulsations. 

Non-radial pulsations in RR Lyr stars were first discov- 
ered by Olech et al. (1999a) in three RRc variables in the 
globular cluster M55. This discovery was quickly followed 
by Olech et al. (1999b) who found one non-radial pulsating 
RRc star in M5, by Kovacs et al. (2000) who reported a few 
such RRc stars in the MACHO data, and by Moskalik (2000) 
who discovered six non-radial RRc and RRab pulsators in 
the OGLE data. 

A theoretical investigation of non-radial pulsations in 
RR Lyr stars was presented by Dziembowski and Cassisi 
(1999). They showed that the largest growth rates are seen 
for radial modes and strongly trapped unstable modes with 
large values of I. The latter do not explain our observations 
because it is difficult to detect pulsations with such large 
values of I (I = 7— 10) in photometric data. The third group 
of modes with quite large growth rates are non-radial pul- 
sations with I — 1 laying in the vicinity of the unstable 
1st overtone mode (see Fig. 2 of Dziembowski and Cassisi, 
1999). The most unstable I — 1 modes always have high 
frequencies (short periods) which is in good agreement with 
our observations. In NGC 6362 there is only one monope- 
riodic RRc star with a period shorter than the periods of 
the RRc stars with non-radial pulsations. The remaining 13 
monoperiodic RRc stars have longer periods. A similar sit- 
uation was observed by Olech et al. (1999a) in M55. Among 



nine RRc stars there is only one star with a period shorter 
than the period of the three non-radially pulsating variables. 

Dziembowski and Cassisi (1999) showed also that the 
relative frequency separation A/// is largest for the I = 1 
modes in the vicinity of the 1st overtone pulsations. For 
the shortest periods A/// is around 0.01 which is in good 
agreement with our observations which give 0.016, 0.017 and 
0.031 for V37, V6 and V10, respectively. 

The B — V color range for RRc variables in NGC 6362 is 
between 0.21 and 0.33. The colors of V37, V6 and V10 are 
0.21, 0.26 and 0.26, respectively, placing these stars near 
the blue edge of the instability strip. The same result was 
obtained by Olech at al (1999a) who found three non-radial 
pulsating RRc stars in M55 with B — V colors between 0.31 
and 0.32. The B — V colors of the RRc stars in M55 are 
between 0.31 and 0.39. 



3.3 Fourier Analysis of the RRab Variables 

Recently Kovacs and Jurcsik (1996, 1997) and Jurcsik 
(1998) have derived formulae that estimate the absolute 
magnitudes, metallicities, intrinsic colors and temperatures 
of RRab stars from the periods, amplitudes and phases as 
measured from a Fourier analysis of the light curves. These 
equations are: 

[Fe/H] = -5.038 - 5.394P, + 1.3450 3 i (7) 
M v = 1.221 - 1.396P - 0.477^i + O.lO30 3 i (8) 
Vo-K = 1.585+1.257P -0.273yli -0.234^31 +0.062041 (9) 
logT e// = 3.9291 - 0.1112(Vo - Ko) - 0.0032[Fe/H] (10) 

where ^41 = <j>4 — 4<f>x (cf. equation 1). 

These equations are valid only for RRab stars with reg- 
ular light curves, i.e. variables with a deviation parame- 
ter Dm smaller than 3 (see Kovacs and Kanbur 1998). In 
our sample, only three variables do not satisfy this condi- 
tion. Table 5 summarizes the results obtained from equa- 
tions (7)-(10). The mean values of the absolute magni- 
tude, metallicity and effective temperature for 14 stars with 
D m < 3 are M v = 0.86 ± 0.01, [Fe/H] = -0.93 ± 0.04 and 
T e ff — 6555 ± 25, respectively. 

In their original paper Jurcsik and Kovacs (1996) in- 
troduced the Dm parameter to distinguish RRab variables 
with regular light curves (D m < 3) from those with irregular 
light curves, i.e. those showing Blazhko behavior (D m > 3). 
In our sample of RRab stars V5, V7, V12 and V13 exhibit 
clear amplitude modulations. Surprisingly, the D m parame- 
ter for these stars is less than 3. This may be a result of the 
large number of observational points in our light curves. The 
Blazhko effect is then averaged over many cycles and the fit 
given in equation (1) is not able to detect the modulation 
since the light curve seems to be regular. 

The estimated mean metallicity [Fe/H] = -0.93±0.04 is 
in very good agreement with the estimates of Zinn and West 
(1984) [Fe/H] = -1.08 and Carretta and Gratton (1997) 
[Fe/H] = —0.96. However, the metallicity computed from 
equation (7) is on the scale of Jurcsik (1995) which is con- 
nected with scale of Zinn and West by the formula: 

[Fe/H] Jufcsik = 1.431[Fe/H] zw + 0.880 (11) 
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Table 5. Parameters for the RRab Variables in NGC 6362 



Star 




P[days] 


A 




Ai 


031 




041 


M v 


[Fe/H] 


T eff 


D m 


VI 





.504800 


15.321 





.376 


5.015 


1 


.412 


0.859 


-1.015 


6572 


1 


.64 








±0.001 


±0 


.001 


±0.009 


±0 


.013 


±0.084 


±0.022 








V2 





.488973 


15.309 


0. 


.422 


5.040 


1 


.357 


0.861 


-0.897 


6637 


1 


.83 








±0.001 


±0. 


.001 


±0.008 


±0 


.012 


±0.084 


±0.022 








V3 





.447110 


15.322 





.449 


5.013 


1 


.293 


0.904 


-0.707 


6726 


2 


.38 








±0.001 


±0. 


.002 


±0.020 


±0 


.038 


±0.082 


±0.035 








V5 





.520730 


15.354 


0. 


.290 


5.220 


1 


.554 


0.899 


-0.825 


6556 


2 


.49 








±0.002 


±0. 


.003 


±0.061 


±0 


.089 


±0.090 


±0.084 








V7 





.521590 


15.323 





.361 


5.113 


1 


.580 


0.852 


-0.974 


6549 


1 


.58 








±0.001 


±0. 


.001 


±0.017 


±0 


.024 


±0.087 


±0.028 








V12 





.532905 


15.278 


0. 


.361 


5.063 


1 


.487 


0.831 


-1.103 


6521 


1 


. 19 








±0.002 


±0. 


.002 


±0.027 


±0 


.040 


±0.086 


±0.040 








V13 





.579987 


15.253 





.362 


5.241 


1 


.864 


0.784 


— 1.118 


6454 


1 


.76 








±0.001 


±0. 


.001 


±0.018 


±0 


.026 


±0.090 


±0.031 








V16 





.525668 


15.346 


0. 


.359 


5.105 


1 


.547 


0.847 


— 1.007 


6541 





.80 








±0.000 


±0. 


.000 


±0.004 


±0 


.006 


±0.087 


±0.018 








V18 





.512876 


15.351 





.344 


5.141 


1 


.541 


0.876 


—0.890 


6571 


1 


.50 








±0.002 


±0. 


.003 


±0.040 


±0 


.055 


±0.087 


±0.056 








V19 





.594506 


15.332 


0. 


.220 


5.522 


2 


.378 


0.860 


—0.818 


6402 


15. 


.26 








±0.000 


±0. 


.000 


±0.010 


±0 


.020 


±0.096 


±0.033 








V20 





.698330 


15.212 


0. 


.163 


5.974 


3 


.459 


0.790 


-0.769 


6227 


1 


.71 








±0.000 


±0. 


.000 


±0.019 


±0 


.058 


0.107 


±0.060 








V25 





.455891 


15.478 





.440 


4.849 


1 


.181 


0.879 


-0.976 


6662 


1 


.51 








±0.000 


±0. 


.000 


±0.004 


±0 


.005 


±0.080 


±0.024 








V26 





.602182 


15.343 


0. 


.227 


5.575 


2 


.507 


0.852 


-0.788 


6395 


2 


.08 








±0.000 


±0. 


.001 


±0.011 


±0 


.024 


±0.097 


±0.035 








V29 





.647743 


15.226 





.184 


5.738 


2 


.314 


0.826 


-0.814 


6366 


15. 


.76 








±0.001 


±0. 


.001 


±0.052 





.114 


0.102 


±0.081 








V30 





.613372 


15.170 


0. 


.316 


5.379 


2 


.075 


0.773 


-1.111 


6396 


2 


.12 








±0.001 


±0. 


.001 


±0.012 


±0 


.018 


±0.094 


±0.029 








V31 





.600230 


15.308 


0. 


.223 


5.617 


2 


.421 


0.861 


-0.720 


6419 


13. 


.89 








±0.001 


±0 


.001 


±0.019 


±0 


.036 


±0.098 


±0.043 








V32 





.497203 


15.330 


0. 


.380 


5.040 


1 


.422 


0.870 


-0.941 


6596 


1 


.50 








±0.001 


±0. 


.001 


±0.012 


±0 


.018 


±0.085 


±0.024 








V34 





.494370 


15.349 


0. 


.263 


5.273 


1 


.658 


0.954 


-0.612 


6599 


2 


.64 








±0.001 


±0 


.002 


±0.031 


±0 


.071 


±0.090 


±0.047 









If we take [Fe/H] = -1.08 from Zinn and West (1984), the 
equation (11) gives [Fe/H] = -0.67 on the Jurcsik (1995) 
scale. The reason for this discrepancy is unknown. Previous 
estimates of [Fe/H] using equation (7) have produced results 
in broad agreement with other estimates (see Kaluzny et al., 
1998, 2000, Olech et al. 1999a). 

3.4 Variable V18 

The light curve of RRab variable V18 shown in Fig. 1 is quite 
noisy. The scatter is observed not only around maximum as 
in Blazhko variables (e.g. V5, V7, V12, V13) but along the 
whole light curve. This may indicate that this variable pul- 
sates in two modes. The non-radial pulsations discovered in 
RRc stars (this paper, Olech et al., 1999a) were also ob- 
served in RRab stars (Moskalik, 2000) and the presence of 
non-radial modes in V18 could explain the scatter in the 
light curve. 

Inspection of the raw light curve suggests an amplitude 
modulation of 0.1 mag to 0.2 mag on a time scale of around 
40 days, perhaps indicating a beating effect. In our opinion 
there is another possible explanation. 

The power spectrum shown in Fig. 8 displays a maxi- 



mum frequency at fo = 1.94979 c/d which corresponds to 
Po = 0.512876 d . There are also many aliases of fo placed at 
fo ± 1, fo + 2 and fo + 3. The AoV spectrum for three har- 
monics also shows a clear variation at fo with weak peaks 
at f ± 1, fo/2, fa/3, fa/4 and (/ + l)/2. 

Both these periodograms are in agreement with a single 
frequency fo but this does not explain the scatter in the light 
curve of V18. One can notice that the doubled period (2 • 
Po = 1.026 d ) differs only 0.026 days from one day and thus 
at each 38 days around midnight we observe only the flat 
part of the light curve which is close to the mean brightness 
of the star. After 20 days around midnight we observe the 
largest change of the amplitude. Thus the modulation of 
the amplitude might not be real but might be caused by 
observing different parts of the light curve at each night. 

The prewhitening of the raw light curve with fo and five 
harmonics gives a noisy light curve with a standard devia- 
tion of 0.07 mag, much larger than the photometric errors. 
The power spectrum of the prewhitened light curve shows 
many peaks with the highest at fo and fo/2 with amplitudes 
of approximately 0.03 mag. 

We conclude that there is no evidence for multiperiodic 
behavior of V18. It is possible that the enhanced scatter 
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detected in the light curve of this star is partially caused by 
the Blazhko effect but this is difficult to check due to the 
proximity of Po to half a day. 

3.5 Reddening and the instability strip of NGC 
6362 

The reddening of RRab variables can be calculated from the 
metallicity, expressed in terms of AS (Preston 1959). From 
Blanco (1992) we have the relation: 

E{B — V) =< B — V >* co .5-o.8) +0.01222AS (12) 
-0.00045(AS) 2 - 0.185P - 0.356 

where < B — V >a>(o.5-o.8) is the observed mean color in 
the 0.5-0.8 phase interval and P is the fundamental period. 
Based on the globular cluster metallicity scale adopted by 
Zinn and West (1984) Suntzeff et al. (1991) derived the fol- 
lowing AS — [Fe/H] relation: 



Table 6. The distance and reddening determinations for NGC 
6362 



[Fe/H] = -0.408 -0.158 AS 



(13) 



Using the equation (13) and adopting [Fe/H] = — 1.0 for 
NGC 6362 we obtain AS = 3.75. 

The average value of E(B — V) calculated in this way for 
our 18 NGC 6362 RRab stars is E(B-V) = 0.08±0.01. This 
value is consistent with recent determination of Piotto et al. 
(1999) who obtained E(B-V) = 0.06±0.03, slightly smaller 
than the value of Alcaino and Liller (1986) who obtained 
E(B-V) =0.10. 

Knowing the colors of RR Lyr variables and the red- 
dening we are able to estimate the boundaries of the insta- 
bility strip in NGC 6362. Our coolest variable is V20 with 
< B > — < V > equal to 0.51 and T e ff computed from 
equation (10) equal to 6227 K. The hottest variable is V14 
with < B > — < V > equal to 0.25 and T e // computed 
from equation (4) equal to 7577 K (we exclude V37 with 
<_B> — < V >— 0.21 due to its noisy light curve). Adopt- 
ing a color excess of E(B — V) = 0.08, we measure the blue 
and red edges of the instability strip in NGC 6362 to be at 
[B - V)o E = 0.17 and (B - V)o E = 0.43. These values are 
consistent with estimates made for other globular clusters 
(Smith 1995). 

3.6 Distance Modulus to NGC 6362 

Piotto et al. (1999) estimate the apparent distance modulus 
to NGC 6362 to be (m - M) v = 14.67 ± 0.20, based on a 
measured mean magnitude of the horizontal branch (HB) 
and the average of a number of calibrations of the relation 
between the absolute magnitude of the zero age horizontal 
branch and the cluster metallicity. 

A comparison of the color magnitude diagram of NGC 
6362 presented in Fig. 2 with that of Piotto et al. (1999) 
shows a small but significant shift in V magnitude. Piotto 
et al. (1999) present a histogram of the distribution in V 
for the HB stars in the color interval 0.5 < B — V < 0.7. 
From this histogram they obtained a mean visual magnitude 
for the HB of Vhb = 15.45 ± 0.03. Fig. 9 shows similar 
histograms (we decided to increase the range of colors to 
0.0 < B — V < 0.85) for all color magnitude diagrams of 
the cluster found in the literature. The panels correspond to 
the work of (from up to bottom) Alcaino and Liller (1986), 



Distance modulus 
(m - M) v 



Color excess 
E(B - V) 



Source 



14.74 

14.67 ±0.20 
14.55 ±0.20 
14.46 ±0.02 
14.59 ±0.03 
14.46 ±0.10 



0.10 

0.06 ± 0.03 



0.08 ±0.01 



Alcaino & Liller (1986) 
Piotto et al. (1999) 
Piotto corrected 

of RRc stars 
log L I Lq for RRc stars 
RRab stars 



Mazur et al (1999), Piotto et al. (1999) and to this work. 
All histograms, except that of Piotto et al., show that the 
mean V magnitude of the HB is equal to 15.33. The 0.12 
mag difference between the Piotto et al. estimate and the 
other data sets in Fig. 9 suggests that the value of apparent 
distance modulus to the NGC 6362 from the Piotto et al. 
work should be (m — M)y = 14.55 ± 0.20 instead of (m — 
M) v = 14.67 ±0.20. 

Knowing the values of Ao and My" for RRc stars (see 
Table 3) we can compute the apparent distance modulus to 
the NGC 6362. It is equal to (m - M)v = 14.46 ± 0.02, 
quite consistent with our corrected value of Piotto's deter- 
mination. On the other hand one should remember that the 
equation for My" is calibrated using Baade-Wesselink lu- 
minosities, which give faint absolute magnitudes for RR Lyr 
stars. We can also use the values of log L/Lq of RRc stars 
to compute the distance modulus. The absolute magnitudes 
My (presented in the last column of Table 3) were calcu- 
lated assuming a value of 4.79 for Mboi of the Sun and using 
a bolometric correction BC = 0.06[Fe/H] + 0.06, adopted 
from Sandage and Cacciari (1990). The resulting distance 
modulus is equal to (m — M)y — 14.59 ± 0.03, in excellent 
agreement with the value presented above. This agreement 
suggests that the theoretical models constructed by Simon 
and Clement (1993) produce the correct values of luminosi- 
ties of RRc stars. 

A similar estimate can be made for the RRab stars. 
Using the data given in Table 5 (14 RRab stars with D m < 
3) we obtain (m — M)v — 14.46±0.10. It is exactly the same 
value as we obtained using Kovacs' expression for My° of 
RRc stars (see section 3.1). The agreement is expected since 
both equations (6) and (8) were calibrated using Baade- 
Wesselink luminosities. 

Table 6 summarizes the distance and reddening deter- 
minations for NGC 6362. 

The values of the distance modulus based on the mean 
magnitude of the HB and that measured from log L/Lq of 
RRc stars are very consistent. It is known that theoretical 
models of RR Lyr stars favor the long distance scale and 
thus large values of distance modulus. On the opposite side 
stand the statistical parallaxes and Baade-Wesselink lumi- 
nosities of RR Lyrae stars which give faint values of absolute 
magnitudes of RR Lyr stars and thus values of the distance 
modulus which are smaller by 0.2-0.3 mag. 

The final conclusion of this section is that the distance 
modulus to NGC 6362 is around 14.45 mag if we prefer the 
short distance scale and around 14.6 if we prefer the long 
distance scale. 
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4 CONCLUSIONS 

We have presented a comprehensive study of 35 RR Lyr 
variables in the globular cluster NGC 6362. Our l/-band 
light curves contain over 1000 measurements and our B- 
band light curves around 260 measurements. The periods of 
the cluster RRc stars are between 0.2462 and 0.3815 days 
with a mean value of 0.295 days. The periods of the RRab 
variables are between 0.4471 and 0.6983 days with a mean 
value of 0.538 days. 

The period-amplitude (log 10 P — Av) diagram for NGC 
6362 shows that RRab variables from this cluster do not fit 
the log 10 P — Av relation for M3 another Oosterhoff type I 
globular cluster, suggesting that the V amplitude for a given 
period is a function of metallicity rather than Oosterhoff 
type alone. 

Fourier decomposition of the light curves was used 
to estimate the basic properties of these variables. From 
the analysis of the RRc stars we obtain a mean mass of 
M = 0.531±0.014 M®, luminosity log L/L® = 1.656±0.006, 
effective temperature T e fj — 7429 ± 20, and helium abun- 
dance Y = 0.292 ± 0.002. The mean values of the ab- 
solute magnitude, metallicity (Jurcsik's scale) and effec- 
tive temperature for 14 RRab stars with "regular" light 
curves are: My = 0.86 ± 0.01, [Fe/H] = -0.93 ± 0.04 and 
T e ff = 6555 ± 25, respectively. 

Using the B — V colors, periods and metallicities of the 
RRab stars, we estimate the value of the color excess for 
NGC 6362 to be equal to E(B - V) = 0.08 ± 0.01. From 
this color excess we derive the unreddened colors of the blue 
and red edges of the instability strip in NGC 6362 to be 
(B - V) BE = 0.17 and (B - V)$ E = 0.43. 

The apparent distance modulus of NGC 6362 calculated 
from the sample of 14 RRab stars is (m—M)v = 14.46±0.10. 
An analysis of the RRc stars based on Kovacs' expression 
for leads to a value of (m - M) v = 14.46 ± 0.02. This 

agreement is expected since both derivations rely on Baade- 
Wesselink luminosities, giving faint magnitudes for RR Lyr 
stars and a short distance scale. On the other hand, a larger 
value of the distance modulus is determined from estimates 
of log L/Lq of the RRc stars - (m - M)v = 14.59 ± 0.03. 
This result is also expected because theoretical models of RR 
Lyr stars produce large values of distance moduli and prefer 
a long distance scale. These latter distance determinations 
are consistent with the value of (m — M)y = 14.55 ± 0.20 
derived from the mean magnitude of the zero age horizontal 
branch. 

Three of the RRc variables presented in Fig. 1 (V37, 
V6 and V10) show a clear modulation of their light curves. 
The power spectra of the light curves of these stars show 
the highest peaks for the 1st overtone radial pulsations. In 
the case of V37 we detect another peak in the vicinity of the 
main frequency which cannot be connected with radial pul- 
sations. In the case of V6 and V10 we detect two non-radial 
modes in the vicinity of the main frequency. Our conclusion 
is that these variables belong to the newly discovered group 
of non-radial pulsating RR Lyr stars (Olech et al., 1999a, 
1999b, Kovacs et al., 2000, Moskalik, 2000). We have also 
noticed that the non-radial pulsation RRc stars lay near the 
blue edge of the instability strip, consistent with the results 
of Olech et al. (1999a) for RRc stars in the globular cluster 
M55. 
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Fig. 1 Light curves of RR Lyr variable stars found in NGC 
6362. The stars are plotted by increasing period. The upper 
light curve is in theV-band, and the lower in the _B-band. 
The B magnitude is reduced by 0.2 mag for clarity. 
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10 A. Olech et al. 

Fig. 1 Continued. 
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11 

Fig. 2 Color-magnitude diagram of NGC 6362. The open 
triangles denote RRab stars and open circles correspond to 
RRc variables. 
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Fig. 3 The phased light curve of variable V25. The observa- 
tional points are plotted with grey open circles. The size of 
these circles is 10 times larger than the photometric errors. 
The solid black line corresponds to the fit computed from 
equation (1). 
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Fig. 4 Period-amplitude diagram for RRab stars with 
D m < 3 (solid triangles), RRab stars with D m > 3 (open 
triangles) and RRc stars (open circles). The solid line rep- 
resents a linear fit to RRab variables in M3 (Kaluzny et al. 
1998) 
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Fig. 5 The power spectrum of the real (top), prewhitened 
(middle) and two times prewhitened (bottom) light curve 
of variable V37. The prewhitened light curves phased with 
periods Pi and P2 are shown on the right side. 
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Fig. 6 The power spectrum of the real (top), prewhitened 
(upper middle), two times prewhitened (lower middle) and 
three times prewhitened (bottom) light curve of variable V6. 
The prewhitened light curves phased with periods Pi, P2 
and P3 are shown on the right side. 
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Fig. 7 The power spectrum of the real (top), prewhitened 
(upper middle), two times prewhitened (lower middle) and 
three times prewhitened (bottom) light curve of variable 
V10. The prewhitened light curves phased with periods Pi, 
P2 and P3 are shown on the right side. 
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16 A. Olech et al. 

Fig. 8 The power spectrum of RRab variable V18. 
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Fig. 9 Histograms of the distribution in V of the HB stars 
in the color interval 0.0 < B — V < 0.85 based on color 
magnitude diagrams obtained by Alcaino and Liller (1986), 
Mazur et al. (1999), Piotto et al. (1999) and in this work. 
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